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Abstract 

We construct analytically the gravity duals of three-dimensional, super Yang- 
Mills-type theories with J\f = 1 supersymmetry coupled to Nf quark flavors. The 
backreaction of the quarks on the color degrees of freedom is included, and corre¬ 
sponds on the gravity side to the backreaction of Nf D6-branes on the background 
of N D2-branes. The D6-branes are smeared over the compact part of the geome¬ 
try, which must be a six-dimensional nearly Kahler manifold in order to preserve 
supersymmetry. For massless quarks, the solutions flow in the IR to an AdS± 
fixed point dual to a Chern-Simons-matter theory. For light quarks the theories 
exhibit quasi-conformal dynamics (walking) at energy scales m q <£< XNf/N, 
with A = < 7 2 m -/V the’t Hooft coupling. 
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1 Introduction 

Three-dimensional, Af = 8 super Yang-Mills (SYM) theory with gauge group SU(N) was 
among the first detailed examples of holographic duality [1]. In three dimensions the gauge 
coupling 5 ym h as dimensions of energy and thus the theory possesses a non-trivial renormal¬ 
ization group (RG) flow. Since the theory is asymptotically free, it provides an example in 
which the ultraviolet (UV) is nicely behaved. The holographic description of this theory is 
obtained from the supergravity solution sourced by N D2-branes in flat space. The amount 
of supersymmetry can be easily reduced to J\T = 1 by replacing the flat space transverse to 
the D2-branes, which is a cone over S 6 , by a GVholononry cone whose base must be a six- 
dimensional nearly Kahler (NK) manifold different from the six-sphere [2], In this case 
the dual gauge theory is presumably a quiver gauge theory. 

Our goal in this paper will be to understand the effect of adding Nr flavors of fundamental 
matter to the theories above. In all cases the flavored theory will be J\f = 1 supersymmetric. 
We will loosely refer to this matter as ‘quarks’ despite the fact that it will include both bosonic 
and fermionic degrees of freedom. We will see that their inclusion leads to interesting infrared 
(IR) dynamics such as the appearance of a Chern-Simons (CS) matter theory in the case of 
massless quarks or of ‘walking’ (quasi-conformal) dynamics in the case of light quarks. 

The addition of quarks on the gauge theory side corresponds to the addition of D6-branes 
on the gravity side [3]. We will work with unquenched quarks, meaning that their backreaction 
on the adjoint (color) degrees of freedom is included. On the gravity side we will therefore 
include the backreaction of the D6-branes on the D2-brane geometry. 

The case when the internal manifold is S' 6 and all the D6-branes are overlapping was 
studied in [4, 5]. The supersymmetric solutions constructed in these references depend on two 
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radial coordinates: the radius transverse to the D2-branes but parallel to the D6-branes, and 
the radius transverse to both the D6- and the D2-branes. Despite their beauty, these solutions 
are difficult to generalize because the dependence on two radii translates into the requirement 
to solve non-linear partial differential equations. For example, it has not been possible to find 
the generalizations corresponding to the introduction of a non-zero temperature or a non-zero 
charge density into the system. 1 * 

In order to reduce the supergravity equations to ordinary differential equations, we will 
smear the D6-branes over the internal manifold [8] (see [9] for a review of this approach). 
This, together with the reduction to first-order equations implied by supersymmetry, will 
allow us to construct the solution essentially analytically. 

The presence of the quarks leaves the UV properties of the gauge theory unmodified, but 
it changes the IR dynamics. In the case of massless quarks, the solution flows to an AdS^ 
fixed point in the IR. We will argue that, in general, the dual gauge description is a CS-matter 
theory. In the particular case in which the NK manifold is CP 3 we will be able to identify 
this theory as the flavored version [10] of the Ooguri-Park solution [11], which is itself an 
AT = 1 deformation of the ABJM [12] theory. In the case of quarks that are light compared 
to the scale XNf/N set by the’t Hooft coupling A = the theory exhibits ‘walking’ or 

quasi-conformal dynamics in the energy range m q <F< XNf/N. 

2 Flavorless solutions and nearly Kahler manifolds 

Since three-dimensional SYM-type theories are asymptotically free, in the UV a perturbative 
description is possible. At an energy scale E ~ A the theories becomes strongly coupled. If N 
is large, then this regime can be described holographically by the gravitational background 
sourced by N D2-branes; in terms of the string coupling and the string length, the Yang-Mills 
coupling is 

£ym=W4- (2-1) 

If the space transverse to the D2-branes is flat space then maximal, J\f = 8 supersymmetry 
is preserved (in three-dimensional language) and the gauge theory is SU(N ) SYM. In order 
to decrease the amount of supersymmetry — as the addition of flavor will do anyhow — the 
stack of branes has to be positioned at the tip of a Ricci-flat cone with reduced holonomy 
[2]. This replaces the transverse flat space, which in polar coordinates is a cone over the 
six-sphere S e . In order to preserve J\f = 1 supersymmetry the cone must have G 2 holonomy. 
The base of a GVcone is a nearly-Kahler manifold. Since properties of NK manifolds will 
play an important role, we now proceed to review them. 

NK geometries appear naturally in the classification of almost Hermitian manifolds [14], 
A NK manifold is an almost Hermitian manifold whose fundamental form J satisfies the 
condition 

3 VJ = dJ, (2.2) 

with V the derivative associated to the Levi-Civita connection. In this terminology, the 
Kahler condition is VJ = 0. Powerful splitting theorems reduce the study of these geome¬ 
tries to six dimensions [15], and we will restrict ourselves to this case in the following. The 
importance of these geometries in physics (and in particular in string-theory) stems from the 
fact that they admit Killing spinors. Indeed, in six dimensions, a manifold admits a Killing 

1 See however [6] for a perturbative, finite-temperature solution near the core of the D6-branes, and [7] for 

a computation of the meson spectrum in the same near-core limit. 
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spinor if and only if it is NK [16]. As a consequence, such a manifold is Einstein with positive 
curvature. 

This property makes NK manifolds suitable geometries to support supergravity solutions 
preserving some amount of supersymmetry. In fact, an equivalent way of defining a NK 
manifold that makes this more apparent is the following. The Riemannian cone 

ds 2 (Cj) = dr 2 + r 2 ds 2 (Ms) (2.3) 

has G 2 holonomy if and only if Ms is NK [17]. As is well known, Go holonomy is the condition 
to preserve minimal supersymmetry in four dimensions starting from eleven-dimensional M- 
theory, in analogy with the SU( 3) holonomy of Calabi-Yau manifolds that is required when 
starting from ten-dimensional string theory. In particular, models with chiral fermions can 
be obtained considering M-theory in the presence of GVcone singularities like (2.3), as shown 
in [18, 19]. 

These internal geometries also appeared in the quest for stabilizing moduli in string com- 
pactifications to Minkowski (and related) vacua. In [20] it was found that massive type 
IIA supergravity admits M = 1, AdS 4 x Ms solutions with Ms being NK. In the context 
of flux compactifications, the most suitable language to analyze the backgrounds is that of 
G-structures. This leads us to another way of characterizing NK geometries, i.e. as six¬ 
dimensional manifolds admitting an SU (3) structure with only Wi non-vanishing among the 
torsion classes in the codification of [21], Possessing SU(3 ) structure implies the existence 
of a globally defined, real two-form J (associated to the almost complex structure) together 
with a globally defined, complex three-form II, satisfying 

J Afl = 0, 

ijAJAJ = ^n An. (2.4) 

From the class of SU (3) structure manifolds, nearly-Kahlerness is selected by the requirements 

dJ = ^ Im [Wi n] , 

dn = WiJAJ, (2.5) 

that is, the torsion classes W 2 , • • •, W 5 corresponding to other SU( 3) representations that 
could occur in (2.5) all vanish. For the purposes of this paper, equations (2.4) and (2.5) will 
be the defining properties of the internal geometries. 

Considerably less attention has been paid to NK manifolds in the context of the gauge/string 
duality, despite the early observation in [2] that they emerge as the transverse space to D2- 
branes preserving J\f = 1 supersymmetry and generalizing S' 6 . In this respect, a parallel can 
be traced with D3-branes in type IIB supergravity and Sasaki-Einstein (SE) manifolds. Just 
like a stack of D3-branes at the tip of a cone over a SE manifold is dual to a gauge theory 
with minimal supersymmetry in four dimensions, a stack of D2-branes at the tip of a GVcone 
as in (2.3) is dual to a gauge theory with minimal supersymmetry in three dimensions. De¬ 
spite this resemblance, there are also some important differences. First, the D2-brane near 
horizon geometry is not AdS, as opposed to the D3-brane case. This means that one cannot 
use the familiar tools of conformal field theories (CFT) to analyze the dual gauge theory. 
And second, compared with Sasakian geometry, NK manifolds are still not well understood 
from the mathematical viewpoint, though several important results are known (for reviews 
see [22, 23]). 
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M d 

C(M d ) 

Gauge theory dual 

G-structure 

Globally defined forms 

NK 

GVcone 

D = 3, N = 1 SYM 

SU{ 3) 

J( 2 ) i ^(3) 

SE 

Calabi-Yau 

D = 4, N = 1 SCFT 

SU{2) 

h(l) 5 J( 2 )) ^( 2 ) 


Table 1: Some key properties of NK manifolds in comparison with SE manifolds. J^) always 
denotes a real two-form. is a complex n-form. r is a real one-form. 


The comparison between NK manifolds in six dimensions and SE manifolds in five dimen¬ 
sions is summarized in the following table: 

The complete list of known, regular, compact, six-dimensional NK manifolds is as follows : 2 


S 6 


g 2 

51/(3) 


CP 3 ~ 


Sp( 2) 

SU( 2) x U{ 1) 


S 3 X S 3 


SU( 2) x SU{ 2) x SU( 2) 
SU{ 2 ) 


F(l, 2) ~ 


SU( 3) 

U( 1) x U( 1) 


( 2 . 6 ) 


All these manifolds are homogeneous. Moreover, any six-dimensional, homogeneous, NK 
manifold is isometric to one of them [27]. We emphasize that the metric on the CP 3 that is 
compatible with the NK structure is not the perhaps-more-familiar Fubini-Study metric used 
in the ABJM construction [12], which is instead Kahler. Regarding CP 3 as an S 2 fibration 
over S' 4 , the NK metric is squashed with respect to the Kahler one. As a consequence, the 
isometry is reduced to Sp( 2) ~ 50(5) C SU( 4). In a similar manner, the NK metric on 
S 3 x S 3 is not the product of the round metrics, so again the isometry is reduced. For more 
details about these cosets we refer the reader to [18]. 


3 Adding flavor 

3.1 Generalities 

We have argued that a stack of N D2-branes placed at the tip of the GVcone (2.3) provides 
the holographic dual — at least in some energy range — to a three-dimensional gauge theory 
with minimal supersymmetry, that is, two real supercharges. When the internal NK geometry 
is taken to be the six-sphere, there is an enhancement of supersymmetry to J\f = 8 and the 
resulting duality was examined in [1]. Still largely unexplored are the detailed held theory 
duals for the remaining cosets in the list (2.6), as well as for other putative NK manifolds to 
be found . 3 

Setting aside what is the specific SYM dual to these gravitational solutions, we will focus 
instead on the inclusion of fundamental matter, to which we will also refer as ‘flavor’ or 
‘quarks’. This is achieved by adding a new set of branes to the gravitational system [3]. 

2 Besides, there are infinite families of examples with orbifold or conical singularities [24, 25, 26]. 

3 See [28] for a proposal for the quiver theory dual to the NK metric on CP 3 . 
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Since we wish to preserve supersymmetry and we want the fundamental degrees of freedom 
to propagate along the 2+1 gauge theory directions (see e.g. the discussion in Section 5.5 of 
[29]), we will introduce an additional stack of Nf D6-branes. The relative orientation between 
the color and the flavor branes is indicated by the following array: 



x 1 

9 

x r 

NK 

D2 

X 

X 


D6 

X 

XXX 

X X 


with each D6-brane wrapping a three-dimensional submanifold inside the NK space transverse 
to the D2-branes. 

The case when all the D6-branes wrap the same submanifold and the internal manifold 
is an S 6 was studied in [4, 5] and generalized to other manifolds in [30]. As explained in 
the Introduction, the supersymmetric solutions constructed in those references solve partial 
differential equations and are therefore difficult to generalize. This technical difficulty can be 
overcome by smearing the flavor branes appropriately [8]. Whenever there is a large number 
of branes, it is possible to distribute them along their transverse directions — the dots in the 
row of the D6-branes in the array above — in such a way that the dependence on the second 
radius disappears. In practice one has traded the partial differential equations for ordinary 
ones. 

In our case we will distribute the D6-branes in a way that preserves the SU (3) structure of 
the NK internal geometry. We will see that this requirement essentially fixes the distribution 
of D6-branes uniquely, and moreover that it is compatible with preservation of supersymmetry. 

The D6-branes will contribute to the energy-momentum tensor and thus modify the metric 
originally sourced only by the D2-branes. The parameter that controls the relative influence 
of fundamental matter with respect to the initial color branes is [9] 


N f 2 

The effective dimensionless coupling g e g is defined as in [1] 

SeS = 91 m NU- 1 = XU- 1 


(3.1) 


(3.2) 


where U is a radial coordinate on the gravity side that is dual to an energy scale in the gauge 
theory. The so-called ‘probe approximation’ on the gravity side, in which the backreaction of 
the flavor branes into the geometry is neglected, is justified when Nf g^g/N is small. On the 
field theory side this corresponds to a ‘quenched’ approximation in which the flavor degrees 
of freedom are treated as probes of the gluon-plus-adjoint-matter-dominated dynamics. If 
Nf 9ef[/N is not small, then the backreaction of the flavor branes on the geometry must be 
included. In the gauge theory this corresponds to a situation with ‘unquenched’ degrees 
of freedom in the fundamental representation. Note that, in the context of the large-lV 
expansion, a necessary condition for the fundamental matter to be unquenched is that the 
limit N —» oo is taken in such a way that Nf/N ^ 0. This way of taking the large- N limit is 
usually referred to as the Veneziano limit. 

It is clear from (3.1) and (3.2) that, for fixed A and non-vanishing Nf/N, there is an 
energy scale, or equivalently a radial position in the geometry, at which Nfg/^/N becomes 
of order unity. Below this scale the backreaction of the D6-branes must be included and we 
expect the geometry to be significantly modified. 
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In the opposite limit, for high enough energies, the backreaction of the flavor branes 
decreases. We therefore expect that the solution of the color branes, dual to pure super 
Yang-Mills, will be recovered in the UV, which is in agreement with the theory being su- 
perrenormalizable. As we will see, all these expectations will be confirmed by our explicit 
solutions. 


3.2 Ansatz 

We will now write down the ansatz for the supergravity fields in our solution and derive the 
corresponding BPS equations, which we will solve in subsequent sections. As indicated by 
the D2/D6 array above, the ten-dimensional geometry consists of three Minkowski directions, 
a radial coordinate and a six-dimensional internal manifold, which is assumed to possess an 
SU (3) structure. This implies that there exist a real two-form J and a complex three-form 
defining the structure and subject to the compatibility conditions (2.4). Given an appropriate 
set of vielbeins e a , a = 1 ,..., 6, at least locally we can write 


J = e 12 + e 34 + e 56 5 

= (e 1 + i e 2 ) A (e 3 + * e 4 ) A (e 5 + * e 6 ) . 


(3.3) 


If in addition the manifold is NK, as we assume, the differential relations (2.5) are verified. 
Without loss of generality, the torsion class can be taken to be real and adjusted to W\ = 2, 
so we have the differential conditions 


dJ = 3ImH, dRefl = 2 J A J, (3.4) 

together with the Hodge duals with respect to the internal metric 

*6 J = —J A J , *6^ = — *12. (3-5) 

The vast simplification implied by these assumptions follows from the fact that we have a 
natural set of forms to employ in the ansatz for the supergravity fields as well as the possibility 
to preserve supersymmetry. 

Let us begin by writing down the ansatz for the supergravity forms in the Ramond-Ramond 
(RR) sector. Since the solution is sourced by D2- and D6-branes, we expect the RR six- and 
two-form field strengths to be non-zero. As explained in Appendix A, the distribution of 
backreacting D6-branes gives rise to a violation of the Bianchi identity for the RR two-form, 

d-Fh = —2 k 2 Td6 H , (3.6) 

where H is the so-called smearing form that indicates how the flavor branes are distributed 
in the internal directions. Intuitively, this is simple to understand. The D6-branes couple 
minimally to the RR potential CV, and hence their presence leads to a source on the right- 
hand side of the equation of motion for its field strength, i.e. we have d * Fg ~ S, where E 
measures the local density and the orientation of D6-branes at any given point. Since by 
definition F 2 = — * F%, this leads to eqn. (3.6). 

Given the forms at our disposal, the simplest ansatz for the RR fluxes reads 

F 2 = Q f J, (3.7) 

F e = —JAJAJ. (3.8) 
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As we will see below, for massless quarks the fact that J is not closed is crucial for the 
consistency of this ansatz with the violation of the Bianchi identity. As usual, the six-form 
sourced by the D2-branes is proportional to the volume form of the internal space. The 
parameters Q c and Qf have dimensions of (length ) 5 and (length ) 1 and are related to the 
number of D2- and D 6 -branes, respectively, or equivalently to the rank of the gauge group 
and the number of flavors, through the quantization condition 


2 ^T Vr / Fa -" ~ (2nl,Y-Vg, J ~ N " ' 

Using (3.7) and (3.8) this immediately gives 

^ ( 2 vr £ s fg s AT ( 2 tt £ s )g s AT 

Qc = - 77 - N > Qf = — i 77 - N f ■ 


(3.9) 


(3.10) 


where the dimensionless quantities Vq and V 2 are the volume of the internal manifold and 
f J, respectively. The computation of V 2 requires the knowledge of an explicit realization of 
the NK structure. We emphasize that, in general, Qf = Qf(r ) can be a function of the radial 
coordinate, as will be the case for massive quarks. Note that, up to numerical coefficients, we 
have that 

Nf 

Qc~A^, Qf-X^l]. (3.11) 

As always in the context of the gauge/string correspondence, the powers of £ s will cancel out 
in the computation of gauge theory observables. 

It is instructive at this point to come back to the analogy with the more familiar SE 
geometries transverse to D3-branes. In this case the flavored solution was found in [31] and 
contains an additional stack of D7-branes, whose backreaction induces a violation of the 
Bianchi identity for the RR one-form: 


d.Fi = —2 k 2 Tb7 S . 


(3.12) 


As shown in Table 1, in a SE geometry we have at our disposal a globally defined, non-closed 
one-fornr. The natural ansatz for the one-form field strength used in [31] was thus 

F 1 = Q f r ] . (3.13) 

Equations (3.12) and (3.13) are the SE analog of our NK equations (3.6) and (3.7). 

Equations (3.4), (3.6) and (3.7) determine the distribution of the flavor branes, i.e. 5, up 
to an arbitrary function Qf(r) as 

“ = ~ 2 K 2 r D6 drAj + 3 Qf Im > ( 3 - 14 ) 

where the ' denotes differentiation with respect to the radial coordinate r. We will see that 
for massless quarks Qf is just a constant, meaning that in this case the branes are smeared 
purely along the directions spanned by Imfh 

We now turn to the Neveau-Schwarz (NS) sector of the theory. Working in string frame, 
we choose to parameterize it as 

ds 2 = h~ 2 dxf 2 + h 1 * e 2x ds 2 (C 7 ) , 


ds 2 (C 7 ) 


e = 


dr 2 + r 2 dsg , 
h* e 3x , 


H = 0. 


(3.15) 
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The internal metric dsg = e a 5 a b e b is that of a NK manifold, associated to the vielbeins 
introduced in (3.3) and normalized so that its Ricci tensor satisfies R a b = 5 5 a b- The G 2 - 
cone metric (2.3) thus appears explicitly in our ansatz. The entire system is specified by two 
functions of the radial coordinate, h(r) and %(r), corresponding roughly to the dilaton and 
the breathing mode (the volume modulus of the internal manifold). Even without allowing 
for deformations of the internal NK manifold, one may have expected the most general ansatz 
preserving the required symmetries to contain an additional function, but this can always 
be integrated in terms of the other two — see Appendix A for details. The virtue of the 
parameterization (3.15) is that the flavorless solution is straightforwardly recovered by setting 
Nf = x = 0 and taking h to be the usual harmonic function of the D2-branes, h = Q c /br 5 . 

At first sight, it may seem overly restrictive to not allow for any deformation of the internal 
metric, specially taking into account that some form of squashing is present in every solution 
with smeared flavor branes discussed in e.g. [9]. At the technical level, this difference is due 
to the fact that in our case there is no natural way of writing a generic NK manifold as a 
fiber over some base. In contrast, every SE manifold is a 17(1) fiber over a Kahler-Einstein 
base, and the backreaction of the flavor branes produces a squashing between the fiber and 
the base. 

One immediate consequence of this difference is that the metrics of our flavored solutions 
possess the same isometries as the unflavored ones. Moreover, when the NK manifold is 
not the six-sphere, the almost complex structure J, and therefore the fluxes in our ansatz, 
preserve all the isometries of the metric [27]. It thus follows that the gauge theories with 
fundamental matter dual to our solutions will be invariant under the same global symmetries 
as the corresponding unflavored ones, the exception being the maximally supersymmetric case 
dual to S' 6 . In this case the metric enjoys a full SO(7) invariance but this is broken down to 
G 2 C SO( 7) by the two-form flux (3.7). This G 2 subgroup can easily be understood as the 
numerator of the coset S 6 ~ G 2 /SU(3). 4 

That this is the preserved symmetry follows for instance from the observation in [32] that 
the forms of the six-sphere NK structure, J and fl, coincide with the G 2 left-invariant forms 
on the coset. The fact that the fluxes preserve a smaller amount of symmetry than the metric 
is consistent with the equations of motion, because the fluxes enter the stress tensor that 
sources the metric only quadratically, and the square of the fluxes is more symmetric than 
the fluxes themselves. 

3.3 BPS equations 

Since we are seeking supersymmetric solutions, our next task is to write down the correspond¬ 
ing BPS equations that the functions in our ansatz must obey. Typically this would require 
studying the fermionic variations of the supergravity fields, but in our case we can bypass this 
by making use of some results in the literature, together with the mathematical properties 
of the background. Indeed, the ingredients needed to write the BPS equations are contained 
in [33]. This reference studies the construction of four-dimensional domain walls in type 
IIA supergravity, understood as solutions possessing Poincare symmetry in three dimensions. 
Minimal supersymmetry is also imposed. NK manifolds are particularly simple examples 
of internal geometries potentially preserving J\f = 1, so they are considered in detail. Our 
ansatz verifies all the assumptions in [33] except for the violation of the Bianchi identity for 

4 This G 2 subgroup should not to be confused with the holonomy of the cone (2.3). The cone over a six¬ 
dimensional NK manifold always has G 2 holonomy, but the isometry group of the NK manifold may be different. 
For example, as mentioned above for CP 3 viewed as a NK manifold the isometry group is Sp( 2) ~ SO(5). 



the two-form. Adapting the results of 
equations in our case read 


ti = -9l e -*X-¥k f ?x h . (3.16) 

Since each of the D6-branes in our solution preserves supersymmetry, we expect that each 
of them wraps a calibrated four-cycle inside the GVcone transverse to the D2-branes. The fact 
that this cone possesses GVholonomy and not just a GVstructure (i.e. that the GVstructure 
is torsion-free, as implied by the fact that the base is NK) means that the corresponding 
associative three-form T and co-associative four-form T = * 7 \k are both closed, and in fact 
they are both calibration forms. In terms of the NK forms and the radial coordinate on the 
cone, the co-associative form is given by 

\k = r 3 dr A Re f! + ^ r 4 J A J. (3-17) 

We thus expect that the full world volume of a given D6-brane obeys a generalized calibration 
condition of the form [34] 

d ( e "* K) = *F 2 = F 8 , (3.18) 

with /C a generalized calibration form proportional to T. We have checked by direct calculation 
that, upon using the BPS equations, this condition is indeed obeyed by 

K, = h* e 4x d 3 x A T . (3.19) 

By definition of a calibration, it then follows that the world volume action of a single D6-brane 
can be written as 

S D6 = -T D6 J {e~*JC-C 7 ) , (3.20) 

where dG 7 = Fg and pull-backs onto the brane’s worlvolume are understood. The general 
results in [35] imply that, in the presence of callibrated branes, the sourced-modified equations 
of motion are implied once the (violated) Bianchi identities are imposed. We have checked 
that this is the case for all the equations of motion listed in Appendix A. 


[33] to account for this fact, we conclude that the BPS 


Qf 


e 2x 
2 ’ 


4 Massless quarks and an infrared fixed point 


We begin by considering the solution for massless quarks which, as we will justify in the 
next section, corresponds to setting Qf to a constant. The system (3.16) of first order BPS 
equations is easily integrable. The general solution, in terms of two integration constants c x 
and Ch, is 


e 


-2x 


c x + 


2 Qf 


h 


( 2 Qf + c \ r ) 


Q c 

315 Q 5 f 


(35 Q) - 20 Q)c x r + 12 Q 2 f c 2 r 2 


+ 



8 Qc \ 

315 Q 5 f J 


c 9 r 9 \ 2 
2 Qf + c x r ) 


8 QfC x r 3 + 8c^ r 4 ) 


(4.1) 
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Without loss of generality we will set c x = 1, since this can be achieved through the rescalings 


c 3 x ^, r —» c x 1 r . 


(4.2) 


The interpretation of Ch is clarified by examining the solution as Qf —>• 0. This limit is smooth 
and results in the unflavored solution corresponding to D2-branes at the tip of a GVcone, 
with x = 0. The dilaton then reads 


= hJ 



(4.3) 


We see that Ch is precisely the constant in the harmonic function of the D2-branes. Setting 
Ch = 1 selects asymptotically flat boundary conditions. Here we will instead set Ch = 0, since 
this implements the decoupling limit that yields the equivalence between the gauge theory 
and the gravity descriptions. In conclusion, after fixing the integration constants, the only 
parameters specifying the solution for massless quarks are the dimensionful’t Hooft coupling 
A and the dimensionless numbers of D2- and D 6 -branes N and Nf, respectively. 

Let us now examine the UV and the IR limits of the solution. The UV regime of the 
gauge theory corresponds to the region r —» oo, in which the leading-order asymptotic form 
of the metric and the dilaton are 


ds 


2 

s 


e 


$ 



da n,2 + 



(dr 2 + r 2 d sf) , 



(4.4) 


This is exactly the solution for N D2-branes, meaning that the addition of flavor does not 
modify the UV properties of the theory at leading order. The first corrections in e.g. the 
dilaton behave as: 


e 


$ 



59 Qf 14657 (QA 2 o (0l 
24 r 2688 \ r ) V r 


(4.5) 


The IR regime of the gauge theory corresponds to the region r —>• 0, in which the metric 
and the dilaton take the form 


ds s = J 2 dx b Jr ^2 d P 2 + \ L ‘ 2 ds 6 > 

* = = 

2 Vs\Q 5 f J 9s2V3\vJ [NjJ ’ 


(4.6) 


where we have changed coordinates via r 3 = Qf p 2 and we have defined 

l = J_ (9A l = Mi (YiY (IV 

3\/3 \Qf) 3\/3 \VJ \Nf) 


(4.7) 


Remarkably, the IR geometry is AdS 4 x NK with radius L. This falls in the general class of 
supersymmetric AdS 4 solutions with SU (3)-structure internal manifold found in [36]. The 
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requirement that type IIA supergravity provides a reliable description, namely the conditions 
that L l s and g s e® -C 1, translate into 


NfCN^Nj. 


(4.8) 


Note that these inequalities require both N and Nj to be large. 

There exists a technical difficulty with uplifting to M-theory type IIA solutions in which 
the Bianchi identity for F 2 is violated, but presumably this uplift can be obtained along 
the lines proposed in [37]. Assuming the usual relation between the radii in ten and eleven 
dimensions and the dilaton, the AdS solution would correspond in M-theory to another AdS 4 
geometry with radius 

L~(NN f )h p , (4.9) 

with l p the eleven-dimensional Planck length. This is always large in the large-A limit, 
thus extending the range of validity of the solution to arbitrary Nf — see a closely related 
discussion in [4], 

We conclude that the solution (4.1) (with c x = 1 and Ch = 0) describes a set of M = 1 su¬ 
persymmetric RG flows (one for each possible internal NK manifold) from a three-dimensional 
SYM theory in the UV to an interacting fixed point in the IR, driven by the addition to the 
theory of massless quarks. We show in Appendix A that the fixed point is approached along 
irrelevant directions corresponding to operators of dimensions A = 6, 11/3. 

Several observations can be made about the gauge theories dual to the solutions above. 
First of all there is the fact that the flow drives the theory to an interacting IR fixed point. 
The existence of non-trivial conformal theories at the IR of a three-dimensional gauge theory 
in the presence of a large number of flavors was first observed in [38]. This was proven in 
an expansion in 1 /Nf to all orders. Here we find that this feature is already present in the 
Veneziano limit, at least for the J\f = 1 cases at hand. To elaborate on this point, we note 
that the running gauge coupling g 2 (g) can be computed introducing a D2-brane probe in the 
flavored background at a fixed radial position r, expanding the D2-brane action to quadratic 
order in the Born-Infeld field strength and matching with the canonically normalized YM 
term. The result for the dimensionless effective coupling is 


9eS = 


g 2 {d)N 


A / A /v,Ayf 

A V 2 N g) 


(4.10) 


To arrive at this equation we have used (2.1) and the fact that the energy scale in the gauge 
theory is related to the radial position in the bulk through 27 t/t = r/l 2 . In the UV and in the 
IR this behaves as 

3 

UV: g 2 eS ~ £, IR: ffeff ~ (j^) 2 ^ • (4-11) 

The UV behavior is that expected of SYM, i.e. (3.2). In between these two asymptotic 
behaviors the coupling attains a unique maximum at 

^ CFT = V 2 ~N ’ (4 ' 12) 

where the subscript is a reminder that the physics below this scale is approximately conformal. 
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At this scale 5 


(4.14) 


2 _ V2 N 
9eS ~ 3 V3N f - 

and the YM /3-function vanishes. This may seem to suggest that the theory reaches a fixed 
point at a finite energy scale. However, this is not true because the IR YM interactions are 
subdominant with respect to the Chern-Simons (CS) interactions generated along the flow, 
which is consistent with the fact that the effective YM coupling actually goes to zero in the 
deep IR. At the perturbative level, the addition of quarks is known to induce a CS term with 
level proportional to the number of flavors running in the loop [40, 41]. At strong coupling, 
this effect can be seen in the non-vanishing of the Wess-Zumino (WZ) action for a probe 
D4-brane that fills out the Minkowski directions and is suitably oriented along the internal 
directions in the flavored background: 


S\vz ~ 


J A AF A 



A A F , 


(4.15) 


where F is the field strength of the gauge field A living on the brane. We see that the presence 
of an F 2 flux induces a CS level proportional to Nf. 

It has been conjectured that the free energy of a three-dimensional (Euclidean) field theory 
placed on the three-sphere is a genuine measure of the number of degrees of freedom and 
verifies appropriate monotonicity theorems [42]. This free energy thus plays a role analogous 
to that of the c-function for two-dimensional field theories or the a-coefficient in the Weyl 
anomaly in four dimensions. For a theory at a fixed point with a holographic description, 
the free energy is proportional to the AdS 4 radius measured in units of the four-dimensional 
effective Planck length. For the IR fixed point (4.6) the free energy scales with the number 
of flavors and the rank of the gauge group as 


F(S 3 ) ~ 


L 2 e~ 2 * 

2k\ 


L 8 e- 2$ 

2k? 


(N 3 N f y , 


(4.16) 


where the dilaton factors account for the fact that the metric (4.6) is written in string frame. 
We see that the ^dependence (the adjoint contribution) matches the peculiar power of CS- 
rnatter theories of the ABJM type. 

The consistent picture that emerges from this analysis is that, below the scale (4.12), the 
CS term dominates over the YM term and the IR dynamics of our models is governed by a 
CS-matter theory. In the case in which the internal NK manifold is CP 3 this can be made 
more explicit because of the connection with the solution by Ooguri and Park (OP) [11], 
which is an JV = 1 deformation of the ABJM [12] theory. Indeed, Conde and Ramallo (CR) 
[ 10 ] constructed solutions corresponding to the addition of backreacting flavor to both the 
ABJM solution and the OP solution. Regarding CP 3 as an S 2 fibration over an S ' 4 base — see 
the paragraph below (2.6) — CR showed that the backreaction produces a relative squashing 
between the fiber and the base, as well as a deformation of the RR form F 2 to accommodate 
the violation of the Bianchi identity. CR parametrize this squashing and this deformation 
with two dimensionless quantities q and 77 , respectively. Given 77 , q = q(rf) is determined by 
a second-order algebraic equation, which results in two branches of solutions. Moving along 
a given branch corresponds to changing the number of flavors. Crucially, the ABJM and the 

s This can be compared with the one-loop result, robust in the large Nf limit (see e.g. the discussion in [39]) 

Ad^gls = — <?eff + Po fleff (4-13) 

for some positive /3q. This vanishes at a non-zero coupling gl B = 1 //3o - 
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OP solutions lie on different branches. By explicit construction of the NK structure on CP 3 , 
it can be shown that our IR fixed point corresponds to a solution with q = 2,rj = — 2 in the 
language of CR, and this solution happens to be on the same branch as the unflavored OP 
solution. In fact, the results of CR show that this solution is obtained by adding Nf = 4|fc| 
flavors to the OP solution. It is interesting to note that the metric on the CP 3 manifold of the 
unflavored OP solution is not even Einstein, yet under the addition of an appropriate number 
of flavors the metric becomes not just Einstein but NK. 

5 Massive quarks and quasi-conformal dynamics 

In the case of massive quarks the D 6 -branes are separated from the D2-branes by a finite 
distance proportional to the quark mass. This translates into the fact that the D 6 -branes 
extend from infinity down to a non-zero minimal value of the radial coordinate r m . By 
Gauss' law this means that the D 6 -brane charge Qf(r ) vanishes for r < r m . In the gauge 
theory this corresponds to the statement that in the IR, i.e. at energies below the quark mass, 
the quarks decouple from the dymanics. In contrast, in the UV the quarks can be treated as 
effectively massless, so Qf(r ) must approach the value that it would have had in the case of 
exactly massless quarks. It is therefore convenient to write 

Qf(r) = QfP( r ), (5-1) 

where Qf is now a constant and p(r) is a dimensionless function that vanishes for r < r m and 
that approaches 1 as r -> oo. In between these two limits p(r ) is monotonically increasing. 

It is important to note that p(r) is not determined dynamically: any function with the 
properties mentioned above is an admissible choice. On the gauge theory side, this freedom 
corresponds to the freedom of adding quarks of many different masses to the theory. The 
supergravity equations can be solved for any p(r ) in terms of an integral for y and a double 
integral for h as 

/ °0 A v 

p ( z ) T 2 > 

h(r) = e~ 3x (c h + Q c . (5.2) 

Taking the massless limit p —>• 1, one can see that the integration constants have the same 
meaning as before, so we can set c x = 1 and Ch = 0. 

There is a family of functions p(r), labelled by r m , such that each member of the family 
describes a theory in which all quarks have exactly the same mass proportional to r m . In 
order to find this family some details about the internal manifold are needed. This means 
that to proceed further we cannot treat all the NK geometries simultaneously. For illustrative 
purposes, we will therefore focus on the six-sphere henceforth. 

The special class of functions corresponding to quarks of equal masses is determined by 
the following consistency condition. By assumption, the full supergravity solution includes 
the backreaction of a large number of flavor branes smeared over the internal manifold by the 
action of a symmetry group — G *2 in the case of D 6 -branes on the S 6 . Each of these branes 
is embedded non-trivially on the internal manifold, but all these embeddings are related by a 
symmetry. This means that the action for the entire set of branes should equal Nf times the 
action of a single one. The mathematical expression of this statement is 

T m [ e~*JC AS = N f T m [ V [JC] m , (5.3) 

J 10 h 


13 



where V denotes the pullback to the brane. The action on the left-hand side is integrated over 
the entire spacetime and depends on the D 6 -brane embedding indirectly only through p(r). 
The action on the right-hand side is only integrated over the seven dimensional submanifold 
occupied by a fiducial D 6 -brane and depends explicitly on the embedding of the brane in 
spacetime. Thus this equation will relate p{r) to the embedding of the brane and hence to 
the quark mass. Note that we should in principle consider the full action instead of just the 
Dirac-Born-Infeld (DBI) part, but this is not necessary in this case because the analysis for 
the Wess-Zumino (WZ) part of the action follows from the DBI analysis by supersymmetry. 

Using our calibration and smearing forms, the left-hand side of (5.3) can be easily com¬ 
puted with the result 

T D6 e~® 1C AH = j eX? ’ 3 (p + ^p'^j d 3 xdr, (5.4) 

where we have explicitly performed the integration over the internal manifold. 

In order to preserve supersymmetry, the three-cycle wrapped by a D 6 -brane inside a NK 
manifold must be calibrated by ReD [43, 44]. In the case of the six-sphere an equatorial 
S 3 C S 6 provides an example of such a cycle, so it is convenient to write the metric on the 
six-sphere as 

dDg = dO 2 + sin 2 0 df2§ + cos 2 6 dSl|, (5.5) 

where ddenotes the metric of the unit-radius n-sphere. The D 6 -brane embedding can then 
be specified as 0 = 0(r). To compute the pullback of the calibration form one needs to use 
the following two results, which are easily obtained using the explicit realization of the NK 
structure presented in Appendix B: 

V [d 3 x A dr A Re D] Dg = sin 4 6 d 3 x A dr A e( 3 ) , 

V [d 3 x A J A J] Dg = 2 cos 6 sin 3 9 9' d 3 x A dr A e( 3 ) , (5.6) 

with prime denoting the radial derivative and e( 3 \ being the volume form of the equatorial 
S 3 . The right-hand side of (5.3) then takes the form 

Toe J e _$ T’[/C ] D6 = T m V 3 J e x r 3 (sin 4 6 > + ^ (sin 4 #)') d 3 xdr, (5.7) 

where V 3 is the volume of the S 3 wrapped by the brane. Comparing (5.4) and (5.7) and using 
the quantization condition (3.10), together with the fact that 12Vg = V 2 V 3 , we arrive at the 
conclusion that 

p(r) = sin 4 0(r). (5-8) 

As anticipated, this equation relates the D 6 -brane embedding Q(r) with the supergravity 
charge density function p(r). The quark mass enters as a boundary condition on the D 6 - 
brane embedding, so our next task is to determine the equation of motion and the boundary 
conditions for the latter. 

The brane embedding is determined by the consistency condition that it solves the equa¬ 
tion of motion of a D 6 -brane probe in the background generated by the D 6 -branes themselves. 
The DBI part of the action for the probe takes the form 

Sdbi = -Tm j = -T m V 3 J e x r 3 sin 3 6 (l + r 2 O' 2 ) * d 3 xdr . (5.9) 
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The WZ part of the action depends on the background RR seven-form. This is determined 
by the condition that dCV = F% = *F- 2 , which is solved by 6 


CV = — Qj e 3x 7 — d 3 x A dr A Re hi. 


(5.10) 


The WZ part of the action is now easily calculated with the result 

Swz = Td6 I V [CV] D6 = — 7d6 V3 I Q/e 3x — sin 4 6 d 3 i'dr . 


(5.11) 


In principle, we should now solve the second-order equations of motion that follow from vary¬ 
ing Sdbi + «S\yz wifi 1 respect to 0 (r). However, in this case we can omit this step because 
we know that the D6-brane embedding is supersymmetric, which means that the brane em¬ 
bedding will be determined by a first-order BPS equation. Moreover, because of the no-force 
condition between different D6-branes, we expect this equation to be the same as for a D6- 
brane probe in an unflavoured background sourced by D2-branes alone. This equation is 
known to be 


9' = 


cot 9 


(5.12) 


and it is easy to verify that the second-order equations of motion are automatically satisfied 
provided (5.12) holds. The solution for the brane embedding is thus 


cos 9 = 



(5.13) 


where r m is an integration constant that determines both the asymptotic behavior of the 
brane embedding and the lowest value of the radial coordinate attained by the brane. From 
the asymptotic behavior we read off that the relation between r m and the bare quark mass 
that enters the UV gauge theory Lagrangian is [45, 46] 


m q 


r m 

2vr i 2 s ' 


(5.14) 


Having solved for the brane embedding we can now use (5.8) to obtain the charge distri¬ 
bution 



(5.15) 


where 0 is the Heaviside theta function. Substituting into (5.2) we obtain the supergravity 
solution for massive quarks: 


=- 2 x _ 


= 


1 1 16 Qj_ 

w 15 rm 


1 + 2 ? (i-IW + Ht) 


2 lr rn \2 . 1 


h = < 


1 + 


16 Qf\ Q, 


15 r m I 5 r 


e- 3x Qcf r °° H 


if r < r n 


if r > r n 


if r <r m 
if r > r m 


(5.16) 


6 Gauge invariance implies that CV is only defined up to exact terms. In view of the calibration condition 
(3.18), we could have also chosen CV = K,. The choice we adopted differs from this one by an exact piece. 
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The constants of integration have been adjusted so that in the UV the D2-brane solution 
is recovered, as for the massless case. In the IR, below the scale set by r m , the D2-brane 
solution is also recovered, but with a different normalization that translates into a finite 
renormalization of the coupling: 

A»=(l + ^ A uv . (5.17) 

V io v m J 


This difference in normalizations is needed to ensure the continuity of the solution at r m . 
The quark mass is similarly renormalized. This can be seen by computing the effective quark 
mass at the decoupling scale, which can be read off from the action of a string stretching 
between the lowest point on a D6-brane and the D2-branes: 


m 


IR 

<7 


1 

2tt£ 2 s 



QttQrr dr 



e x dr 


r m f-\ \ 16 Qf\ 2 

27 t £2 15 r m ) 


(5.18) 


Note that to)) 1 < m q . 

The solution exhibits two qualitatively different regimes depending on the value of the 
ratio Qf/r m , which in terms of gauge theory parameters may be written as 


Qf A Nf fj, c ft 
r m m q N m q 


(5.19) 


with /r c ft the scale introduced in (4.12). If m q -C /r CFT then the theory first reaches the region 
at ~ fj ,cft in which the physics is described by a conformal CS-matter theory, as described 
in Section 4, and only at a much lower scale m q it ‘realizes’ that the quark mass is non-zero. 
Thus in this case the theory exhibits ‘walking’ at scales m q <C ^ <C ir CFT , i.e. the physics is 
approximately conformal in this window. In contrast, if m q > /r CFT , then the quarks decouple 
from the dynamics before their presence can drive the theory to an approximately conformal 
phase and the walking region disappears. These two regimes are clearly seen in the behavior 
of the dilaton displayed in Fig. 1, obtained by numerically integrating (5.16). The outcome 
only depends on the ratio between scales, which are clearly visible in the plot. The lower one 
is approximately m q , that we keep fixed. The’t Hooft coupling A is also common to all curves 
in order to have the same UV asymptotics. The different curves are only distinguished by the 
value of Nf/N , which increases from top to bottom. In the bottom curves a quasi-conformal 
region in which the dilaton is approximately constant is clearly visible. As we decrease Nf/N 
and move to the top curves, the size of the walking region decreases until it disappears all 
together. 


6 Outlook 

We have constructed analytic solutions dual to three dimensional SYM theories with dynam¬ 
ical flavors of arbitrary mass. The simplicity of the solutions, in particular the fact that they 
can be found by solving ordinary differential equations, opens the door to several interesting 
extensions including the addition of temperature and/or a quark density. In the first case we 
expect that, as the ratio m q /T is varied, the solutions will exhibit the backreacted version of 
the phase transitions between Minkowski and black hole embeddings uncovered in [47, 48, 49]. 
In these second case we can make contact [50] with the results of [51] in the limit in which 
the quarks mass becomes sufficiently large. 

7 In fact the solution at r m is not just continuous but C 2 . 
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Figure 1: Dilaton as a function of the energy scale for the solutions with massive quarks. The values 
of r m and g s on the gravity side, or equivalently of m q and A on the gauge theory side, are the same for 
all curves. Instead, the ratio Nf/N increases from top to bottom. Specifically, we take Q c = and 
Qf = r m 10" with n = {—2,0, 2,4, 6} from top to bottom, which translates into Nf/N = 10 n m q /A. 
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A Conventions 

In this appendix we provide the details needed to reproduce the solutions in the main text 
together with some additional material, like the reduction of the system to a four-dimensional 
action. 


A.l Ten-dimensional action and equations 

We work with type IIA supergravity in string frame, enlarged with adequate sources. The 
total action is: 

S = SllA T ^sources = S'nS + <SrR T S'sources . (A.l) 

The first part is the action of type IIA supergravity. The Neveu-Schwarz sector is 


£ns = 


l 

2 ^ 


e 


- 2 $ 


* 1 + 4d4> A *d<F 


\h A *h) , 


(A.2) 


with H = d B solving the Bianchi d H = 0. The RR piece, containing the kinetic terms for 
the two- and four-forms plus a topological interaction (whose explicit form will not be needed 
here), reads 


Srr 


1 

2k 1 


— -Fb A *F2 — -F 4 A *F 4 ) + S t0 p 


(A.3) 
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In the absence of sources, the field strengths are 

F 2 = dCi, F 4 = dC 3 - H A Ci , (A.4) 

solving the Bianchi identities 

dF 2 = 0 , d F 4 = H A F 2 . (A.5) 

On the other hand, the action for the sources is the sum of the smeared DBI and WZ terms 8 

Sources = -To6 J (e“* A - C 7 ) A H , (A.6) 

where JC is the calibration form, given essentially by the induced metric on the brane, and Ci 
is the Hodge dual of Ci, with the convention 


F 2 = — * Fg =>■ Fg = *F 2 

F 4 = *Fq => F 6 = -*F 4 . 


(A.7) 


This new coupling between the D 6 -branes and CY is translated into a modification of the 
Bianchi for F 2 (which coincides with the equation of motion for Fg) that now reads 

d F 2 = —2k 2 Td6 E • (A. 8 ) 

The equations of motion for the lower forms are unmodified by the sources 

d (e" 2$ *H)-F 2 A *F 4 - ^F 4 A F a = 0 

d * F 2 + H A *F 4 = 0 (A.9) 

d * F 4 + H A F 4 = 0. 

In contrast, the equation of motion for the dilaton receives an additional contribution from 
the DBI action 

R * 1 + 4d * d4> — 4d4> A *d<f> — A *H — k 2 Td 6 e # /C A S = 0 . (A.10) 

Einstein’s equations are also modified to reflect the presence of the explicit sources 

Rmn + 2V M Viv$ - \h mrs H n rs = T% + Tj^ rces , (A. 11) 

where is the usual energy-momentum tensor of the RR sector 9 


rTiIIA 

1 MN 


mr{F 2 )n R + 


t^{F4)mr 1 r 2 r 3 (F 4 ) N RiR2R3 


8 Here and in all subsequent equations we assume that the gauge invariant combination T = 2iva'dA + B 
vanishes, where A is the Born-Infeld field on the branes. 

9 For arbitrary n-forms uj n and we have defined the contraction symbol as 1 ' n - 
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The contribution from the sources is 

WW** = ^ [smn Ho (*/C) - S .MRS (*IC) N RS ] . (A.13) 

In deriving the last equation we used the equation of motion for the dilaton, as well as the 
facts that 1C A 3 = — [3j (*/C)] * 1 depends on the metric only through the vielbeins and that 

5 _ 1 A 6 / \ 

~ 2 ■ (A ' 14j 

The WZ term in the action does not contribute to this equation, since it is topological, its 
only effect being the modification of the Bianchi identity for the two-form. 


A.2 BPS equations and reduction to four dimensions 

The conditions for having J\f = 1 solutions preserving three-dimensional Poincare invariance 
with a NK internal geometry in type IIA were derived in [33]. We can capitalize on their 
results to write the BPS equations for our system simply by allowing for a suitable violation 
of the Bianchi identity for the 2-form. In order to make contact with their notation, here we 
use the following ansatz for the NS sector 10 

e 2A 

ds s = — dx i,2 + P { dz2 + ds l) > 


e* = S- 

T 


H = 0, 


(A.15) 


where dsg is the metric of the internal NK manifold, normalized to R a b = 5 5 a b- The RR-forms 
are taken as in Eq. (3.7). Then, the set of BPS equations that can be gathered from [33] by 
including the effect of the flavor is 


t 

P 


A' 


Q c 

1 

Qf p 

2 

P T ' 

2 r 

Qc 

1 

3 Q f p 

2 

P T 

2 r 

Qc 

1 

3 Qf p 

4 

P T 

4 r 


(A-16) 


As usual, the equation for the warp factor can be integrated in terms of the rest of the 
functions. In this case it can be written as 


. . 3 1 , 

A = A 0 + - z + - log r , 


(A.17) 


where A$ is a content-free integration constant that can be absorbed into a rescaling of 
the Minkowski coordinates. Thus, the information of the system is contained in just two 
functions, which justifies the form of our ansatz in the bulk of the paper. The dictionary with 
the variables used there is 


p = r 2 h 2 e 2x , 


r = r 3 h 2 , 


e = r. 


(A.18) 


10 With the exception that we have changed to a more convenient radial variable pdr t here = rdzhere- 
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It is also convenient to have a reduced action from which the BPS equations follow. The 
reduction of (massive) type IIA supergravity on an arbitrary NK manifold was performed in 
[52], The resulting theory is an JV = 2 gauged supergravity in four dimensions. From the 
full collection of modes respecting the NK structure, our solutions excite a very restricted 
subset comprising the dilaton and the breathing mode, that is, the volume modulus of the 
internal space. On top of that we have to allow for the deformation coming from the flavor 
branes, which was not considered in [52]. The reduced four-dimensional action that captures 
the dynamics of our brane intersection reads 11 * * 


1/*/ 3 2 3 ^ ' 

S' 4 = 7^2 J ( R * 1 - 7^2 dp A *d p - ^2 dr A *dr - V * l - Q' f * 1 


= 2^2 / (-R * 1 — Gij d</>* A *dc/A — V * 1 — Q' f * 1 J , 


(A.19) 


where the potential is 


2 /Ur 4 2 t 4 t 6 pr z 


(A.20) 


The first piece, quadratic in Q c , descends from the kinetic term of the four-form in ten 
dimensions. The last one is due to the curvature of the internal manifold. This two terms 
will appear in the reduction on any arbitrary Einstein manifold, not necessarily NK. The 
additional two terms containing Qf are exclusive of NK geometries in the presence of flavor. 
The one quadratic in the number of flavors is inherited from the kinetic term of the two-form, 
while the linear term is due to the DBI action. Notice that, at least when Q'f = 0, this is a 
consistent truncation, in the sense that any solution to the action (A.19) can be uplifted to a 
type IIA solution in the presence of sources as described in the previous section. 

The potential can be derived from the superpotential 


W = 


p ( 12 r - 3Q f p) - Q c 

Api T 2 


by means of the standard relation 


(A.21) 


V = (d- 2) [(d - 2) G ij diWdjW - (d - 1) W 2 ] . (A.22) 


Here d is the dimension of the bulk, four in our case. Using the domain wall ansatz 

dsl = e 2A dx\ 4 _ 2 + dz 2 (A.23) 

all the equations of motion are solved by the first order system 

A' = ±W (\> V = T(d~ 2) G ij djW. (A.24) 

The superpotential has an extremum at 

T = \ (QcQf ) 5 , (A.25) 

11 The unusual term Q), proportional to the radial derivative of the flavor charge distribution, is needed in 

the massive case. Varying with respect to Qf gives the calibration condition. Here and in the following we 

will use a four-dimensional radial coordinate, related to the ten-dimensional one by rdzio = d« 4 . 


l (Qc 

!> 3 U/ 
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which is of course an extremum of the potential as well. This uplifts to a supersymmetric 
AdSi geometry, in the class found in [36], that governs the IR of the more general solutions 
to the full set (A. 16). From the derivatives of the potential we can compute the masses of the 
scalars around the AdS point. They turn out to be 

m 2 L 2 = 18, y , (A.26) 

corresponding to irrelevant operators of dimensions 

A = 6, y . (A.27) 

B Nearly Kahler structure of S G 

In this appendix we provide the details for constructing a NK structure on the six-sphere, 
needed to obtain the massive quarks solution. Contrary to other NK manifolds, there exist 
infinitely many such structures compatible with the round metric on the sphere [27]. In [44] 
a constructive algorithm is given by considering the cone over S' 6 , that is M 7 , as the space of 
octonions. We will simply state the final result in a convenient set of coordinates. 

Let us write the metric of the six sphere as 

dflg = & 6 2 + sin 2 6 df2§ + cos 2 6 df2|, (B.l) 

where dflo and dfl 2 are respectively the metrics of a two and three sphere 

dfl| = da 2 + sin 2 aida|, dO§ = d/3 2 + sin 2 f3\ (d/?| + sin 2 /32d/3f) . (B.2) 

Take the following functions of the angles: 


u 1 

= cos 6 cos a \, 

u 2 

= cos 9 sin a\ cos «2 , 

u 3 

= cos 9 sin sin 02 , 

u 4 

= sin 9 cos (3 1 , 

u 5 

= sin 9 sin /3\ cos (5 2 , 

u e 

= sin 9 sin j3\ sin cos /3s , 

u 7 

= sin 9 sin /3\ sin @2 sin ^3 , 




dictated by the embedding of the unit-radius sphere into flat space, = 1. In terms 

of the one-forms v n = d u n , the metric can be written as ^ v n v n = dOg, while the almost 
complex structure reads 

J = y}{v 2 3 +V V + V 56 )+U 2 (v 45 +V 67 -V 13 )+U 3 (v 12 +V 57 -V 4(i ) 

+ u 4 (v 36 - v 17 - v 25 ) + u 5 (v 24 - v 16 - v 37 ) + u 6 (v 15 - v 27 - v 34 ) 

+ u 7 (V 4 + v 26 + v 35 ) . (B.4) 

On the other hand the complex three-form has real part 

Re 0 = u 1 (u 257 - u 246 - u 345 - u 367 ) + u 2 (u 146 + u 356 + u 347 - u 157 ) 

+ V 3 (u 145 + u 167 - u 247 - u 256 ) + u 4 (v 237 + v 567 - v 126 - u 135 ) 

+ u 5 {v 127 + v 236 + u 134 - u 467 ) + u 6 (V 24 + u 457 - u 137 - u 235 ) 

+ u 7 {v 13e -v 125 -v 234 -v 356 ) , (B.5) 
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as well as imaginary part 


Im n = V 123 + v 156 + v 147 + v 245 + v 267 - v 346 + v 357 . (B. 6 ) 

As mentioned in the main text, a BPS D 6 -brane has to wrap a three-dimensional submanifold 
X calibrated by Re 17. This implies Imn|v = 0, as shown in [44]. An example is the equatorial 
S 3 C S e given in our coordinates by 6 = 7t/2. 


References 

[1] N. Itzhaki, J. M. Maldacena, J. Sonnenschein and S. Yankielowicz, “Supergravity and 
the large N limit of theories with sixteen supercharges,” Phys. Rev. D 58 (1998) 046004 
[hep-th /9802042]. 

[2] B. S. Acharya, J. M. Figueroa-O’Far rill, C. M. Hull and B. J. Spence, “Branes at conical 
singularities and holography,” Adv. Theor. Math. Phys. 2 (1999) 1249 [hep-th/9808014]. 

[3] A. Karch and E. Katz, “Adding flavor to AdS / CFT,” JHEP 0206 (2002) 043 [hep- 
th/0205236], 

[4] O. Pelc and R. Siebelink, “The D2 - D 6 system and a fibered AdS geometry,” Nucl. 
Phys. B 558, 127 (1999) [hep-th/9902045]. 

[5] S. A. Cherkis and A. Hashimoto, “Supergravity solution of intersecting branes and 
AdS/CFT with flavor,” JHEP 0211 (2002) 036 [hep-th/0210105]. 

[ 6 ] M. Gomez-Reino, S. G. Naculich and H. Schnitzer, “Thermodynamics of the localized 
D2-D6 system,” Nucl. Phys. B 713, 263 (2005) [hep-th/0412015]. 

[7] J. Erdmenger and I. Kirsch, “Mesons in gauge / gravity dual with large number of 
fundamental fields,” JHEP 0412, 025 (2004) [hep-th/0408113]. 

[ 8 ] F. Bigazzi, R. Casero, A. L. Cotrone, E. Kiritsis and A. Paredes, “Non-critical holog¬ 
raphy and four-dimensional CFT’s with fundamentals,” JHEP 0510 (2005) 012 [hep- 
th/0505140], 

[9] C. Nunez, A. Paredes and A. V. Ramallo, “Unquenched Flavor in the Gauge/Gravity 
Correspondence,” Adv. High Energy Phys. 2010, 196714 (2010) [arXiv: 1002.1088 [hep- 
th]]. 

[10] E. Conde and A. V. Ramallo, “On the gravity dual of Chern-Simons-matter theories 
with unquenched flavor,” JHEP 1107 (2011) 099 [arXiv: 1105.6045 [hep-th]]. 

[11] H. Ooguri and C. S. Park, “Superconformal Chern-Simons Theories and the Squashed 
Seven Sphere,” JHEP 0811 (2008) 082 [arXiv:0808.0500 [hep-th]]. 

[12] O. Aharony, O. Bergman, D. L. Jafferis and J. Maldacena, “N =6 superconformal Chern- 
Simons-matter theories, M2-branes and their gravity duals,” JHEP 0810 (2008) 091 
[arXiv:0806.1218 [hep-th]]. 

[13] N. Seiberg, “Notes on theories with 16 supercharges,” Nucl. Phys. Proc. Suppl. 67 (1998) 
158 [hep-th/9705117], 

[14] A. Gray and L. M. Hervella, “The sixteen classes of almost Hermitian manifolds and 
their linear invariants,” Ann. Mat. Pura Appl. 123 (1980) 35 

[15] P. A. Nagy, “NK geometry and Riemannian foliations,” Asian J. Math. 6 (2002) 481 
[math.DG/0203038], 


22 


[16] R. Grunewald, “Six-dimensional Riemannian manifolds with a real Killing spinor,” Ann. 
Global Anal. Georn. 8 (1990) 43 

[17] C. Bar, “Real Killing spinors and holonomy,” Comm. Math. Phys. 154 (1993) 509 

[18] M. Atiyah and E. Witten, “M theory dynamics on a manifold of G(2) holonomy,” Adv. 
Theor. Math. Phys. 6 (2003) 1 [hep-th/0107177]. 

[19] B. S. Acharya and E. Witten, “Chiral fermions from manifolds of G(2) holonomy,” hep- 
th/0109152. 

[20] K. Behrndt and M. Cvetic, “General N = 1 supersymmetric flux vacua of (massive) type 
IIA string theory,” Phys. Rev. Lett. 95 (2005) 021601 [hep-th/0403049]. “General N=1 
supersymmetric fluxes in massive type IIA string theory,” Nucl. Phys. B 708 (2005) 45 
[hep-th /0407263]. 

[ 21 ] S. Chiossi and S. Salamon, “The intrinsic torsion of 51/(3) and G 2 structures,” 
[math.DG /0202282]. 

[22] M. Verbitsky, “An intrinsic volume functional on almost complex 6 -manifolds and NK 
geometry,” Pacific J. Math. 235 (2008) 323. 

[23] C. P. Boyer and K. Galicki “Sasakian Geometry, Holonomy, and Supersymmetry,” 
[math.DG/0703231], 

[24] M. Cvetic, G. W. Gibbons, H. Lu and C. N. Pope, “Bianchi IX selfdual Einstein metrics 
and singular G( 2 ) manifolds,” Class. Quant. Grav. 20 (2003) 4239 [hep-th/0206151]. 

[25] K. Behrndt, G. Dall’Agata, D. Lust and S. Mahapatra, “Intersecting six-branes from 
new seven manifolds with G(2) holonomy,” JHEP 0208 (2002) 027 [hep-th/0207117]. 

[26] M. Fernandez, S. Ivanov, V. Munoz and L. Ugarte, “Nearly hypo structures and com¬ 
pact NK 6 -manifolds with conical singularities,” J. London Math. Soc. 78 (2008) 580. 
[math.DG/0602160]. 

[27] J. B. Butruille, “Homogeneous NK manifolds,” [math.DG/0612655]. 

[28] A. Loewy and Y. Oz, “Branes in special holonomy backgrounds,” Phys. Lett. B 537 
(2002) 147 [hep-th/0203092]. 

[29] J. Casalderrey-Solana, H. Liu, D. Mateos, K. Rajagopal and U. A. Wiedemann, 
“Gauge/String Duality, Hot QCD and Heavy Ion Collisions,” CAmbridge University 
Press 2014, [arXiv:1101.0618 [hep-th]]. 

[30] B. S. Acharya, F. Denef, C. Hofman and N. Lambert, “Freund-Rubin revisited,” hep- 
th/0308046. 

[31] F. Benini, F. Canoura, S. Cremonesi, C. Nunez and A. V. Ramallo, “Unquenched flavors 
in the Klebanov-Witten model,” JHEP 0702 (2007) 090 [hep-th/0612118]. 

[32] D. Cassani and A. -K. Kashani-Poor, “Exploiting N=2 in consistent coset reductions of 
type IIA,” Nucl. Phys. B 817 (2009) 25 [arXiv:0901.4251 [hep-th]]. 

[33] M. Haack, D. Lust, L. Martucci and A. Tomasiello, “Domain walls from ten dimensions,” 
JHEP 0910 (2009) 089 [arXiv:0905.1582 [hep-th]]. 

[34] J. Gutowski, G. Papadopoulos and P. K. Townsend, “Supersymmetry and generalized 
calibrations,” Phys. Rev. D 60, 106006 (1999) [hep-th/9905156]. 

[35] P. Koerber and D. Tsimpis, “Supersymmetric sources, integrability and generalized- 
structure compactifications,” JHEP 0708 (2007) 082 [arXiv:0706.1244 [hep-th]]. 

[36] D. Lust and D. Tsimpis, “Supersymmetric AdS(4) compactifications of IIA supergravity,” 
JHEP 0502 (2005) 027 [hep-th/0412250]. 


23 



[37] J. Gaillard and J. Schmude, “The Lift of type IIA supergravity with D6 sources: M- 
theory with torsion,” JHEP 1002, 032 (2010) [arXiv:0908.0305 [hep-th]]. 

[38] T. Appelquist and D. Nash, “Critical Behavior in (2+l)-dimensional QCD,” Phys. Rev. 
Lett. 64 (1990) 721. 

[39] J. Polchinski, “Introduction to Gauge/Gravity Duality,” arXiv:1010.6134 [hep-th], 

[40] A. J. Nierni and G. W. Semenoff, “Axial Anomaly Induced Fermion Fractionization and 
Effective Gauge Theory Actions in Odd Dimensional Space-Times,” Phys. Rev. Lett. 51, 
2077 (1983). 

[41] A. N. Redlich, “Gauge Noninvariance and Parity Violation of Three-Dimensional 
Fermions,” Phys. Rev. Lett. 52 (1984) 18. “Parity Violation and Gauge Noninvariance 
of the Effective Gauge Field Action in Three-Dimensions,” Phys. Rev. D 29 (1984) 2366. 

[42] D. L. Jafferis, I. R. Klebanov, S. S. Pufu and B. R. Safdi, “Towards the F-Theorem: N=2 
Field Theories on the Three-Sphere,” JHEP 1106, 102 (2011) [arXiv:1103.1181 [hep-th]]. 

[43] J. Gutowski, S. Ivanov and G. Papadopoulos, “Deformations of generalized calibrations 
and compact nonKahler manifolds with vanishing first Chern class,” Asian J. Math. 7 
(2003), 39 [math.DG/0205012] 

[44] P. Koerber and L. Martucci, “D-branes on AdS flux compactifications,” JHEP 0801 
(2008) 047 [arXiv:0710.5530 [hep-th]]. 

[45] A. Karch and E. Katz, “Adding flavor to AdS / CFT,” JHEP 0206, 043 (2002) [hep- 
th/0205236], 

[46] M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, “Meson spectroscopy in AdS 
/ CFT with flavor,” JHEP 0307, 049 (2003) [hep-th/0304032], 

[47] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik and I. Kirsch, “Chiral symmetry 
breaking and pions in nonsupersymmetric gauge / gravity duals,” Phys. Rev. D 69, 
066007 (2004) [hep-th/0306018]. 

[48] D. Mateos, R. C. Myers and R. M. Thomson, “Holographic phase transitions with fun¬ 
damental matter,” Phys. Rev. Lett. 97, 091601 (2006) [hep-th/0605046]. 

[49] D. Mateos, R. C. Myers and R. M. Thomson, “Thermodynamics of the brane,” JHEP 
0705, 067 (2007) [hep-th/0701132], 

[50] A. F. Faedo, A. Kundu, D. Mateos, C. Pantelidou and J. Tarrio, “Super Yang-Mills with 
Compressible Quark Matter,” to appear. 

[51] A. F. Faedo, A. Kundu, D. Mateos and J. Tarrio, “(Super)Yang-Mills at Finite Heavy- 
Quark Density,” JHEP 1502, 010 (2015) [arXiv:1410.4466 [hep-th]]. 

[52] A. -K. Kashani-Poor, “Nearly Kaehler Reduction,” JHEP 0711 (2007) 026 

[arXiv:0709.4482 [hep-th]]. 


24 



